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Magnetospheric filter effect for Pc 3 Alfv6n mode waves 

X. Zhang, 1 R. H. Comfort, 1 D. L. Gallagher, 2 J. L. Green, ^ Z. E. Musielak, 1 and 
T. E. Moore 2 

Abstract We present a ray-tracing study of the propagation of Pc 3 Alfvdn mode waves origi- 
nating at the dayside magnetopause. This study reveals interesting features of a magnetospheric 
filter effect for these waves. Pc 3 Alfvdn mode waves cannot penetrate to low Earth altitudes 
unless the wave frequency is below approximately 30 mHz. Configurations of the dispersion 
curves and the refractive index show that the gyroresonance and pseudo-cutoff introduced by the 
heavy ion 0 + block the waves. When the 0 + concentration is removed from the plasma compo- 
sition, the barriers caused by the O'*' no longer exist, and waves with much higher frequencies 
than 30 mHz can penetrate to low altitudes. The result that the 30-mHz or lower frequency 
Alfvdn waves can be guided to low altitudes agrees with ground-based power spectrum observa- 
tions at high latitudes. 
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Introduction 

Studies have shown the magnetopause to be an important 
wave source region for Alfv6n mode Pc 3 waves. Guided prop- 
agation of Alfv6n waves in the magnetosphere has generally 
been discussed in terms of so-called toroidal and transverse 
poloidal modes of oscillation [ Dungey , 1954, 1967]. It has 
been established that many of the disturbances in the Earth’s 
magnetosphere are a response to variations in the solar wind. 
Much evidence has shown that plasmas can penetrate through 
the Earth’s bow shock to ionospheric altitudes in the cusp 
region [Heikkila and Winningham , 1971; Frank , 1982]. There 
are several possible mechanisms for the generation of Alfv6n 
mode waves. One is the perturbation of fluid elements at the 
magnetopause by a sudden change in solar wind pressure. This 
change will move the boundary from its equilibrium position, 
and the disturbance would propagate in the form of Alfv6n 
waves [ Sibeck , 1990]. Another wave generation mechanism is 
the surface wave model, involving the Kelvin-Helmholtz 
instability at the magnetopause [, Miura , 1984, Yumoto et ai , 
1984]. A third wave generation mechanism is the penetration 
by upstream ULF waves of the collisionless bow shock 
region, across the magnetosheath and the magnetopause. 
Studies have shown that Pc 3 pulsations are related to activity 
in the ion foreshock region of the solar wind upstream of the 
Earth’s bow shock and magnetosphere [Odera, 1986; Arnoldy 
et aly 1988; Engebretson et al , 1990]. 

In a previous paper [Zhang et aly 1993] we examined propa- 
gation properties of compressional Pc 3 waves generated in 
the upstream dayside magnetopause using a ray-tracing analy- 
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sis. The results show that the magnetosphere filters the high- 
frequency components of Pc 3 compressional waves, allowing 
the low-frequency component waves to penetrate to low alti- 
tudes. The penetration of Pc 3 compressional waves is, to a 
great extent, controlled by the He + -0 + relative concentration. 
In this paper we present preliminary results of a ray-tracing 
study of Pc 3 Alfvdn mode waves (20 mHz < / < 100 mHz) 
started from the dayside magnetopause. We demonstrate that a 
similar magnetospheric filter effect also applies to Pc 3 
Alfvdn mode waves. We then compare the ray-tracing results 
with ground observations from the South Pole Station and 
McMurdo, Antarctica (at invariant latitudes of -74° and -79°, 
respectively), reported by Engebretson et al. [1989, 1990]. 

Models for This Study 

In a previous paper [Zhang et aly 1993] we introduced our 
ray-tracing program and the plasma and magnetic field models. 
Briefly, the basic ray-tracing equations are from Haselgrove 
[1954], and this ray-tracing program is based on the Stix 
[1962] cold plasma index of refraction. In this study we use a 
modified ray-tracing code from Green et al. [1977], The origi- 
nal ray-tracing code was written by Shawhan [1967]. The 
magnetic field model used in this study is from the empirical 
model of Mead and Fairfield [1975]. The distribution of 
plasmaspheric plasma density is from the empirical density 
model of D, L. Gallagher et al. (Composite empirical model of 
the inner magnetosphere, submitted to Journal of Geophysical 
Researchy 1995), and plasma composition is based on the 
observations of thermal ion distributions by Horwitz et al 
[1986] and Comfort et al [1988]. For this study we use 75% 
H + , 20% He + , 5% 0 + throughout the magnetosphere. 

Characteristics of Dispersion Curves and the 
Refractive Index for Pc 3 Alfven Mode Waves 

Figure 1 shows cold plasma dispersion relation curves for 
the model magnetopause region in the Pc 3 wave frequency 
range. In the dispersion curves we mark two types of wave 
modes, Alfv6n mode and fast mode. Figure 2 shows contours of 
the refractive index for the transverse Alfv6n mode wave in the 
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Figure 1. Cold plasma wave dispersion curves for condi- 
tions at the equatorial day side magnetopause (11 R E geocen- 
tric). 


meridian plane for a frequency of 100 mHz. Comparing this 
topology of the refractive index of Pc 3 Alfv6n mode waves 
with the dispersion relation (Figure 1), we can identify three 
resonance regions in Figure 2 which correspond to the three 
ion cyclotron resonances for H + , He + , and 0 + for Pc 3 Alfv6n 
mode waves. For a fixed frequency, heavier ions will have a 
cyclotron resonance closer to the Earth, at higher magnetic 
field strengths, than lighter ions. The H + resonance therefore 
is found at the largest radial distance from the Earth (in the 
polar region), the 0 + resonance is closest to the Earth, and the 
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He + resonance is ih between. The other two intense gradient 
regions correspond to the plasmapause and to the boundary 
between the closed field and polar regions. For decreasing 
wave frequency the resonances will occur at larger radial dis- 
tances. This is seen in Figure 3 which shows the contours of 
the refractive index for a wave frequency of 40 mHz. Note that 
the 0 + resonance no longer encircles the Earth, but is open at 
low latitudes. The significance of this becomes apparent be- 
low. 

In order to understand the behavior of Alfv6n mode waves, 
we need to discuss some new features besides cutoffs and 
gyroresonances and their corresponding characteristic frequen- 
cies. In the three-ion-species dispersion relation, there are two 
locations in which the fast mode waves and the Alfv6n mode 
waves have the same phase velocity at zero wave normal 
angle. The corresponding frequency is called a cross-over fre- 
quency, defined by Smith and Brice [1964]. For three-ion- 
species plasmas we have two cross-over frequencies, which we 
denote by F CT \ and F ct 2 . These cross-over frequencies occur be- 
tween the H + -He + and He + -0 + resonances, respectively. There 
are two other special locations in the dispersion curves (Figure 
1). The first of them is located in the frequency range between 
the He + -0 + cutoff frequency and F 0 +, the 0 + cyclotron fre- 
quency, and the second between the H + -He + cutoff frequency 
and F He +, the He + cyclotron frequency. At these two special 
locations the dispersion curves show steep slopes, similar to 
resonances and cutoffs. We will refer to them as pseudo-cutoffs 
(or pseudo-resonances), since their properties are similar to 
cutoffs or resonances for waves. Their corresponding frequen- 
cies are called bi-ion-hybrid frequencies [Rauch and Roux , 
1982], which we refer to as F bil and F bi2 , respectively. 

We note in Figure 3 (and by inference from Zhang et al . , 
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Figure 2. Log scale contours of the refractive index for 
Alfv6n mode waves in the noon-midnight meridian plane (100 
mHz). 


Figure 3. Log scale contours of the refractive index for 
Alfv6n mode waves in the noon-midnight meridian plane (40 
mHz). 
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[1993]) that when the wave frequency is lower than about 40 
mHz, the 0 + resonance band becomes open at low latitudes. 
This frequency (40 mHz) is near the “critical frequency” for Pc 
3 Alfv6n mode waves. Similar to that defined by Zhang et al 
[1993], this critical frequency is defined as the frequency for 
which the 0 + resonance location just becomes tangent to the 
magnetopause at the equator. The value of this frequency de- 
pends on the magnitude of the geomagnetic field at the mag- 
netopause. As such, the critical frequency will vary with solar 
wind conditions and may vary with the model of the geomag- 
netic field being employed. As the field intensity increases the 
critical frequency will also increase. 

Waves with frequencies below the critical frequency are 
located leftward of Fq+ in the dispersion relation curves (Figure 
1). Since the general tendency for these waves is to travel 
toward the intense magnetic field strength region along field 
lines, the dispersion curves move to the right relative to the 
wave frequency. Therefore Pc 3 Alfv6n mode waves below the 
critical frequency (or below the local 0 + gyrofrequency) are 
free from all cutoff or resonance (including pseudo-cutoff or 
pseudo-resonance) barriers. These waves are able to penetrate 
to the Earth’s ionosphere. Unlike the situation of cutoffs for 
fast mode waves, the locations of the gyroresonances do not 
vary with the plasma relative concentrations. However, if we 
eliminate an ion species, say 0 + , the corresponding resonance 
disappears. The influence of pseudo-resonances is discussed 
further below. 

Ray-Tracing Examples of Alfvln Mode Waves 

In the Figures below, we present families of ray paths for Pc 
3 transverse AlfvSn mode waves. All rays are launched from 



Figure 4. Meridian ray paths for 100-mHz Alfv6n mode 
waves starting at the equatorial dayside magnetopause (11 R E 
geocentric). 



Figure 5. Meridian ray paths for 40-mHz Alfvdn mode 
waves starting at the equatorial dayside magnetopause (1 1 R E 
geocentric). 

the magnetic equator on the dayside magnetopause. Rays are 
launched in the noon-midnight magnetic meridian plane at 
varying initial wave normal angles. We use a solid line for 
rays with WKB condition [Zhang et al . , 1993] equal to or less 
than 0.1 and a dashed line when the WKB condition is greater 
than 0.1. A solid line means that the ray tracing well satisfies 
the WKB approximation; a dashed line gives a warning for the 
WKB approximation. 

Figures 4, 5, and 6 show rays of transverse Alfv6n waves 
with frequencies of 100 mHz, 40 mHz, and 30 mHz. We display 
rays launched with a range of initial wave normal angles from 
0° to 180°, at 10° intervals. We can see that Pc 3 Alfv6n mode 
waves are guided by the magnetic field line, as is characteristic 
of Alfv6n waves and required by the index of refraction. From 
Figure 4 the 100-mHz rays are guided along magnetic field 
lines to the increased magnetic field strength region so that 
each of the gyrofrequencies increases. In the dispersion rela- 
tion (Figure 1) we can see that the 100-mHz wave frequency is 
between Fq+ and F Hc +. As the magnetic field increases along 
the wave trajectory, the dispersion curves in Figure 1 move to 
the right. At the location where F = F cr2 , the waves undergo a 
polarization reversal and are converted from left-hand 
polarized to right-hand polarized waves. As the waves move 
farther along the dispersion curve, they encounter the pseudo- 
cutoff at the point F = F bi2 , the second bi-ion-hybrid 
frequency, and waves are reflected. If we compare this ray plot 
with the contours of the refractive index for both fast mode 
[Zhang et al , 1993] and Alfv6n mode (Figure 2) waves for this 
same frequency, we find that the location of reflection for 100 
mHz in Figure 4 is between the He + -0 + cutoff and the 0 + 
gyroresonance. This position just corresponds to the location 
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Figure 6. Meridian ray paths for 30-mHz Alfv6n mode 
waves starting at the equatorial dayside magnetopause (11 R E 
geocentric). 

of the pseudo-cutoff in the dispersion relation (Figure 1). In 
this case the range of frequencies at which waves are reflected 
back to the magnetopause depends on the relative concentra- 
tions of He + and 0 + . Reflection will occur at higher frequencies 
as a result of an increasing ratio of 0 + to He + concentrations. 

The situation for 40-mHz rays in Figure 5 is quite different. 
This frequency (40 mHz) is just slightly below F 0 + in the dis- 
persion relation (Figure 1). These AlfvSn mode waves propa- 
gate along an L shell into a region of decreasing magnetic 
field strength (a feature of the realistic dayside magnetic field). 
In the dispersion relation (Figure 1) the gyrofrequencies move 
leftward, showing that the wave will propagate into the 0 + 
gyroresonance and be absorbed. If we overlay the contours of 
the refractive index for 40 mHz for Alfv6n mode waves (Figure 
2) with the ray trajectory traces (Figure 5), we can see that 
these rays stop at the location where they encounter the 0 + 
gyroresonance. 

The 30-mHz rays in Figure 6 are below the critical fre- 
quency. The wave source is in a location below F 0 + in the dis- 
persion relation (Figure 1). When these waves propagate 
along magnetic field lines, they narrowly miss (individual 
rays encounter) the 0 + gyroresonance. After these rays pass 
the region of decreasing magnetic field strength they enter a 
region of increasing magnetic field strength, so that the wave 
locations move leftward along dispersion curves in Figure 1 
(in the dispersion relation (Figure 1) the gyrofrequencies move 
rightward). Therefore when the rays escape from the 0+ 
gyroresonance, they are also free from all the cutoff or reso- 
nance (including pseudo-cutoff or pseudo-resonance) barriers 
at frequencies above Fq+. Therefore most 30-mHz rays are able 
to be guided along magnetic field lines to the Earth’s iono- 
sphere. 


Since an 0 + concentration in the magnetosphere introduces 
cutoff and resonance (including pseudo-cutoff or pseudo- 
resonance) barriers in the magnetosphere in the Pc 3 frequency 
range, it is of interest to check the case where there is no 0 + 
concentration. Figure 7 shows Alfvdn mode rays at 100 mHz 
with no 0 + , where the concentration of H + is 80% and He + is 
20%. Because the 0 + concentration is removed, all the cutoffs 
and resonances (including pseudo-cutoffs or pseudo- 
resonances) at frequencies below F Hc + in the dispersion rela- 
tion (Figure 1) are gone. There are no more barriers for 100- 
mHz Alfvdn mode waves, so they are able to penetrate to the 
Earth’s ionosphere. This plot further supports our suggestion 
that 0 + is the main factor which keeps Pc 3 transverse Alfvdn 
waves from reaching the ionosphere. 

Discussion and Summary 

There have been extensive experiments and theoretical 
reports [e.g., Anderson et al , 1990] on the role which Pc 3 
waves (20-100 mHz) play in energy transport from the solar 
wind to the inner magnetosphere. However, the fundamental 
question as to how wave energy is transported from the solar 
wind in the foreshock magnetosheath to the low-altitude 
magnetosphere remains unclear. Our previous ray-tracing 
study [Zhang et al . , 1993] revealed important features regard- 
ing transport of wave energy from the magnetopause into the 
inner magnetosphere through Pc 3 compressional waves. This 
study reveals magnetospheric wave energy transport through 
the other wave mode branch, Pc 3 Alfv6n mode waves. 

The penetration of upstream ULF magnetic pulsations 
across the magnetopause was studied by Troitskaya et al. 
[1971], Hoppe and Russell [1983], and Yumoto et al. [1984]. 



Figure 7. Meridian ray paths for 100-mHz Alfv6n mode 
waves (in a magnetosphere without 0 + ) starting at the equa- 
torial dayside magnetopause (11 R E geocentric). 
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Yumoto et al. [1984] suggested that the magnetopause wave 
source of high-latitude Pc 3 pulsations can be either surface 
waves excited by solar wind-driven instabilities or upstream 
waves in the Earth’s foreshock. Observations support the idea 
that wave energy could be transported via cusp field lines from 
the magnetopause to the ionosphere through Pc 3 Alfv6n mode 
wave propagation [Bol'shakova and Troitskaya , 1984; 
Plyasova-Bakounina et al. , 1986; Morris and Cole , 1987; 
Engebretson et al ., 1989], Our ray-tracing study of Pc 3 Alfv6n 
mode waves provides a more concrete picture of how wave 
energy at the magnetopause might be transported to the high- 
latitude ionosphere and also what constraints on this transport 
might be imposed by the magnetospheric plasma medium. 

This study shows that the dominant behavior of the Earth’s 
magnetosphere is like that of a low pass filter for transverse 
Alfv6n waves generated at the equatorial magnetopause. Figure 
8 shows schematically the characteristics of the filter action 
of the magnetosphere for Pc 3 transverse Alfv6n mode waves. 
The vertical coordinate is wave frequency, and the horizontal 
coordinate is the penetration altitude. This Figure shows that 
when the wave frequency is below about 30 mHz, Pc 3 Alfv6n 
mode waves are able to penetrate to the \-R E (ionosphere) 
region, so the passband to the ionosphere for Pc 3 transverse 
Alfv6n mode waves is below 30 mHz. There is also a stop band 
in the characteristics of the filter. The stop band is from about 
40 mHz to about 60 mHz, for which there are no Alfv6n mode 
waves. The reason for the existence of the stop band is that 
between the gyrofrequency and bi-ion cutoff, AlfvSn mode 
waves do not exist. At somewhat higher frequencies AlfvSn 
mode waves can propagate into the magnetosphere before 
being turned back by the He + -0 + ion-ion resonance. Just 
above this band and above the lower frequency transmission 
band, the Alfvdn mode wave encounters either the He + or the 
0 + gyroresonance and is absorbed. Of course, the frequency 
spectrum which is actually observed at a given location 
depends on the source spectrum as well as the propagation 
characteristics. 

The results from our ray-tracing study match very well with 
high-latitude ground station observations. Engebretson et al. 
[1989, 1990] reported that Pc 3 wave activity was observed at 
South Pole Station, which is located at a latitude near the 



Magnetopause Ionosphere 

Radial Distance 

Figure 8. Schematic characteristics of the filter action of 
the magnetosphere for Pc 3 Alfv6n mode waves. 


nominal footpoint of the dayside cusp/cleft region, and 
McMurdo, Antarctica (at invariant latitudes of -74° and -79°, 
respectively). During selected days in March and April 1986, 
Pc 3 activity correlated with the low interplanetary magnetic 
field (IMF) cone angle, suggesting that this Pc 3 activity was 
related to upstream wave activity. The power spectra observed 
at these two high-latitude observatories have a steep drop-off 
in spectral density toward high frequencies. The cutoff fre- 
quency in the passband of the power spectrum is near 30 mHz, 
which is consistent with the results from our ray-tracing stud- 
ies. As noted above, the observed drop-off in wave frequency 
could be associated with the source mechanism for the waves. 
To determine whether source characteristics could also be 
responsible would require comparison of wave spectral obser- 
vations at both high and low altitudes along the appropriate 
field lines. However, it would be an interesting coincidence if 
both source and propagation mechanisms gave similar obser- 
vational consequences at low altitudes. 

From the above analysis we can conclude that the cold 
plasma approximation is sufficient to explain a low-pass 
effect for Alfv6n mode waves propagating from the equatorial 
magnetopause to the ionosphere. The passband is located 
below about 30 mHz; this frequency is a function of the 
magnetopause magnetic field strength. That means that 
increased solar wind pressure should raise the magnetic field 
strength there and increase the maximum frequency of Pc 3 
waves observed on the ground. The characteristics of the 
magnetospheric filter are largely determined by the magnetic 
field strength at the magnetopause and along the mag- 
netopause L shell. However, removal of the 0 + concentration 
in our plasma model dramatically enlarges the passband of the 
filter. In that case the magnetosphere lets almost the full 
frequency range of Pc 3 transverse Alfv6n waves pass to the 
ionosphere. 

We anticipate that Alfv6n mode waves will also be able to 
reach the ionosphere for very low, nonzero concentrations of 
0 + . A measure of the likelihood of tunneling between the He + - 
0 + bi-ion cutoff frequency and the 0 + gyrofrequency is taken to 
be when the required tunneling spatial distance is equal to the 
wavelength of the Alfv6n waves. For the situation shown in 
Figure 4, 100-mHz waves are reflected away from the Earth at a 
distance of about 8 R E . At the bi-ion cutoff, Alfv6n waves have 
a wavelength of about 0.1 Rg. We find the tunneling distance 
will be the same as the wavelength when the ratio of the 0 + to 
He + relative concentrations falls below 0.04. If H + has a con- 
centration of 75%, as used in this work, then tunneling is ex- 
pected to become important for concentrations of 0 + less than 
about 1% of the total density. 

Since the cleft ion fountain is known to pump heavy ions 
into the magnetosphere, the question arises as to whether 
heating and associated upwelling O 4, could effectively change 
the Pc 3 wave characteristics observed in the ionosphere. The 
answer is probably not. First, to make a change, the dayside 
cleft flux tubes would have to initially contain no 0 + at high 
altitudes, and then suddenly have 0 + flow up to sufficiently 
high altitudes (low magnetic field intensity) that its gyrofre- 
quency could resonate with the wave. Since ions flowing up 
the cleft flux tubes are typically convecting plasma toward the 
nightside, and the heavier ions are convected more strongly 
than the lighter ions due to the mass spectrometer dispersive 
effect, it is not likely that 0 + would reach the necessary alti- 
tudes on the dayside to have a significant effect, from this type 
of mechanism. 

Of course this mechanism is not limited to the dayside. At 
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any location where the waves are propagating, where 0 + is 
present, and where the magnetic field magnitude is at the right 
level, significant effects on the propagation could be produced 
by this filter mechanism. If it is operating, the observable 
effect would be a broader spectrum of waves (at the high end) at 
magnetospheric altitudes than at ionospheric altitudes. 

The characteristics of the Pc 3 Alfvdn mode waves are quite 
analogous to those of Pc 3 compressional fast mode waves. 
Our previous study [Zhang et al , 1993] also reveals the low- 
pass feature of the magnetosphere to the Pc 3 compressional 
fast mode waves. The mechanism of this low-pass feature for 
Pc 3 compressional fast mode waves is also attributed to the 
existence of the 0 + concentration in the magnetosphere, 
which along with the He + concentration produces the He + -0 + 
cutoff barrier to the propagation of the Pc 3 compressional 
fast mode waves. The difference between the propagation 
characteristics of Pc 3 Alfvdn mode waves and Pc 3 compres- 
sional fast mode waves is that the Alfv6n mode waves could 
propagate to the high-latitude ionosphere along the magnetic 
field line, while the compressional fast mode waves could 
penetrate to the low-latitude plasmasphere, if their frequencies 
were all lower than corresponding critical frequencies. Another 
difference is that the propagation characteristics of Pc 3 
Alfv6n mode waves are less sensitive to the relative concentra- 
tions of helium and oxygen ions than are the Pc 3 compres- 
sional fast mode waves. 
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